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Abstract 

Existing data on the (unpolarised) transverse structure function of the proton are ana- 
lyzed for large values of the Bjorken variable x. The leading- twist and a phenomenolog- 
ical higher-twist contributions are simultaneously determined from a power correction 
analysis of the Nachtmann moments for values of the squared four-momentum trans- 
fer between ~ 1 and 20 (GeV/c) 2 . The results obtained adopting the next-to-leading 
order approximation and those including the effects of soft gluon resummation are com- 
pared. The sensitivity of the extraction of large- x higher twists to high-order radiative 
corrections as well as to the value of a s (M§) is illustrated. 
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The experimental and theoretical investigation of the behaviour of parton distribution 
functions (PDF) for large values of the Bjorken variable x as well as the extraction of higher 
twist terms in the Deep Inelastic Scattering (DIS) structure functions of the nucleon has 
recently received a renewed attention [0-0. The main goals are: i) the determination of 
multiparton correlations in the nucleon 0-0, and ii) the investigation of ratio of d- to 
u-quark distributions in the vicinity of the threshold region x ~ 1 ^ |I| . In this way one may 
test not only models of the structure of the nucleon (see, e.g., [0), but also the predictions of 
perturbative QCD at large values of x, which can be relevant for the experimental programs 
at high-energy facilities, like CEBAF, HERA and future LHC. 

The aim of this letter is to address the issue of the relevance of high-order radiative 
corrections in the extraction of both the leading- and the higher-twist contributions from the 
(unpolarised) DIS structure function of the proton for large values of the Bjorken variable 
x. To this end: i) we perform power correction analyses of the Nachtmann moments of the 
transverse proton structure function, M n (Q 2 ), constructed in |J starting from world data 
in DIS kinematics as well as in the nucleon-resonance production regions for values of the 
squared four-momentum transfer Q 2 between ~ 1 and 20 (GeV/c) 2 and ii) we compare 
the results obtained adopting a fixed-order calculation with those including the effects of 
soft gluon resummation (SGR) || [|. The Nachtmann moments M n (Q 2 ) will be analyzed 
in terms of the following twist expansion (cf. |J) 
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where fi n (Q 2 ) is the twist-2 contribution, while the pQCD evolution of the twist-4 [twist- 



6] is accounted for by the term (a s (Q 2 )/a s (Qo)) 7 " ' [(oc s (Q 2 ) / oc s (Qo))' y " '] an effective 
anomalous dimension 7^ [7^] (cf. also 0) and the parameter affi [a^] represents the 
overall strength of the twist-4 [twist-6] term at the chosen scale Qq. In what follows we will 
limit ourselves to moments with order n > 4, mainly because the second moment is only 
slightly sensitive to the large-x region and it turns out to be almost totally dominated by the 
leading twist (cf. 0). Moreover, the effects of soft gluon resummation for the second moment 
turn out to be quite small (see later on Fig. 1). Finally, since at large x the evolution of the 
singlet quark and gluon operators decouple, we consider the non-singlet (NS) evolution for 
the full twist-2 term (cf., e.g., f|). 

For the fixed-order calculation of the leading twist we consider the next-to-leading 

b Note that existing SLAC data up to x ~ 0.98 as well as the elastic peak contribution are included in the 
construction of the Nachtmann moments performed in ||. Moreover, the inclusion of data in the nucleon- 
resonance production region is a relevant point in view both of the extraction of multiparton correlations 
(target-dependent higher twists) and of parton-hadron duality arguments (see again ||). 

c We want to remind that the Nachtmann moments of Ref. || have been obtained using available inter- 
polation formulae for evaluating the nucleon structure functions at fixed values of Q 2 . In this respect new 
systematic measurements are needed to confirm the Q 2 -behaviour of the Nachtmann moments found in j^] . 
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(NLO) order approximation, viz. (cf. [|I0| 
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where are the one-loop NS anomalous dimensions, 7 i 1 ' iVS) are the two-loop NS anoma- 
lous dimensions, /?o = 11 — 2JV//3, /3i = 102 — 38iV//3 (with iVj representing the number of 
active flavours) and A^ LO) is the n-th moment of PDF at the scale Q 2 . In Eq. (@) C^ ) 
is the NLO part of the quark coefficient function, that in the MS scheme reads as 
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where Cp = (iV^ — 1)/(2N C ) — > 4/3 and Sfc(n) = z3j=i For large n (corresponding to 

the large- x region) the coefficient C^ L °^ is logarithmically divergent; indeed, since S\(n) = 



7^ + log(n) + 0{l/n) and S 2 (n) 
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The physical origin of the logarithm and double logarithm terms in Eq. (|3]) is the 
mismatch among the singularities generated in the quark coefficient function by the virtual 
gluon loops and the real gluon emissions, the latter being suppressed as the elastic peak 
corresponding to the threshold x = 1 is approached. In other words at large x the relevant 
scale is no more Q 2 , but Q 2 (l — x) ]11[ and the usual Altarelli-Parisi evolution equation 
should be accordingly modified [|12[ . The presence of the above-mentioned divergent terms 
at large n would spoil the perturbative nature of the NLO approximation (Q) (as well 
as of any fixed-order calculation) and therefore the effects of soft gluon emissions should 
be considered at all orders in the strong coupling constant a s . To this end one can take 
advantage of resummation techniques, which show that in moment space soft gluon effects 
exponenziate || |13|| . Thus, the moments of the leading twist including the SGR effects, 
acquire the following form (cf. |14], |T5|): 
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the key quantity of the soft gluon resummation and reads as (cf. |L3|, [TJJ 
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where A[a a ] = C F a s /ir + C F K(a s /ir) 2 /2, B[a s ] = -3C F a s /2n with K = C A (67/18-7r 2 /6) - 
10T fl JV>/9, C A = N c -> 3 and Tr = 1/2. Explicitly one has 

G n (Q 2 ) = Hn)G 1 (\ n ) + G 2 (K) + 0[a k s \n k -\n)} (7) 

where X n = /3oa s (Q 2 )\n(n) / 4ir and 
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It is straightforward to check that in the limit A n << 1 one has G n (Q 2 ) — > a s {Q 2 )C\^ L L QQ/2'n, 
so that at NLO Eq. (|) reduces to Eq. (@). 

In || the Nachtmann moments M n (Q 2 ) have been analyzed adopting the NLO ap- 
proximation (0) and using the value a s (M|) = 0.113 derived in |16|]. However, both a recent 
re-analysis of the combined BCDMS and SLAC data |J and the updated PDG sum- 
mary []T7|| of various determinations of a s (M§) clearly indicates a larger value for ct s (M§). 



Therefore, in order to investigate not only the impact of high-order perturbative corrections, 
but also the relevance of the value of a s (M|) we will consider in our analyses two cases: 
a s {Ml) = 0.113 and a s {M 2 z ) = 0.118. 

The results obtained for the ratio of the quark coefficient function calculated within 
the SGR and at NLO, namely 
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are reported in Fig. 1. It can be seen that soft gluon effects are quite small for the second 
and the fourth moments, while for n > 6 they increase significantly as n increases and exhibit 
a remarkably sensitivity to the value of a s (M|), particularly for Q 2 ~ few (GeV/c) 2 . 

Using the twist expansion (Jl]), the Nachtmann moments M n (Q 2 ) of Ref. |J have 
been analyzed adopting the SGR expression ([5]) for the twist-2 moments and the values 
a s (Mf) = 0.118 and Qq = 10 (GeK/c) 2 . Our results, including the uncertainties for each 
twist term generated by the fitting procedure^, are reported in Fig. 2 for n > 4. It can 
be seen that: i) the extracted twist-2 term yields an important contribution in the whole 
Q 2 -range of the present analysis, and ii) the Q 2 -behaviour of the data leaves room for a 

d In our analyses the twist-2 parameter A^ GR ^ [or A^ 1 " ^ in case of NLO analyses] as well as the higher 
twist parameters , , a4 , 7« 6 ^ are simultaneously determined from a x 2_mm i m ization procedure (cf . 
H ) . We have checked that within the errors generated by the fitting procedure the same value for the leading 
twist parameter Aif GR ^ [ J 4^ JVZ/0 ' ) ] is obtained by limiting the range of our analysis only at large Q 2 [namely 
Q 2 > n (GeV/c) 2 ] without including higher twist terms in Eq. (fth. 
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higher twist contribution positive at large Q 2 and negative at Q 2 ~ 1 -j- 2 (GeV/c) 2 ; such a 
change of sign requires in Eq. (JJ) a twist-6 term with a sign opposite to that of the twist-4. 
As already noted in ||, such opposite signs make our total higher twist contribution smaller 
than its individual terms. 

We have then carried out twist analyses of the Nachtmann moments M n (Q 2 ) adopting 
the NLO approximation (0) and using also the value a s (M|) = 0.113. The comparison of 
the results obtained in this way, is shown in Figs. 3-4 for the leading twist fi n (Q 2 ) and for 
the total higher twist contribution WT n (Q 2 ), defined as (cf. Eq. ([J)) 
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From Fig. 3 it can clearly be seen that both soft gluon effects and the value of a s (M§) can 
change drastically the evolution of the leading twist at Q 2 ~ few (GeV/c) 2 |. Correspondingly 
(see Fig. 4), the role of the higher twists is significantly reduced both by the effects of the 
soft gluon resummation and by the value adopted for a s (M§). On the one hand, within 
the NLO approximation, the extracted higher twists for n > 6 can exceed 50% of the 
leading twist at Q 2 ~ few (GeV/c) 2 both for a s (M 2 z ) = 0.113 and a s (M 2 z ) = 0.118. On the 
other hand, within the SGR approach, the total higher twist contribution can be still very 
significant at Q 2 ~ few (GeV/c) 2 if the value a s (M|) = 0.113 is adopted, whereas it does 
not exceed ~ 30% of the leading twist for Q 2 > 2 (GeV/c) 2 , when the updated PDG value 
a s (M z ) = 0.118 is considered. Therefore, we can state that the NLO approximation for the 
leading twist lead to a strong overestimate of the extracted higher twists at large x. 

Before closing, we want to comment briefly on Eq. ( [370] ) as well as on various models 
used in the literature to describe higher twists at large x (see Refs. JL|-|| and fl6fl). As for 
the latter ones, an interesting way to describe the x-shape of higher twists is the renormalon 
approach [[3], in which both 1/Q 2 and 1/Q A power corrections appear as the result of an 
ambiguity in the summation of a specific class of high-order radiative corrections. However, 
the main drawback of the renormalon approach is the fact that the overall strength of the 
renormalon contribution cannot be predicted theoretically and should be fixed phenomeno- 
logically (see Refs. |1|, ||, 0). We can say that renormalons parameterize our ignorance about 
high-order perturbative effects, whereas the approach employed in this work (the soft gluon 
resummation) provides a way to calculate the contribution of a broad class of high-order 
radiative corrections without introducing any additional unknown parameter. 

At present, a rigorous QCD-based treatment of higher twist operators, including their 
mixing under the renormalization-group equations, is not completely available. Therefore, 
one has to resort to phenomeno logical ansatze, like those used in f|, ^, 16 1 and in the present 
work, i.e. Eq. ( |T0"D (cf. also 0, §])• We want to point out that the ansatze used in 0, |6], 116] , 



which by the way are limited only to a twist-4 term, correspond to special cases of Eq. (|1( 
As a matter of fact, the perturbative evolution of the twist-4 is neglected in 0, while it 
is assumed to be the same as the one of the leading twist in [16[]; this corresponds in Eq. 



e The impact of the results shown in Fig. 3 on the large-x behaviour of PDF will be presented elsewhere. 
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(|ICi|) to assume = and = 7^ 5 , respectively. On the contrary, in || as well as 
in the present work we have found non-vanishing twist-4 (and twist-6) effective anomalous 
dimensions, which differ significantly from the values of the twist-2 anomalous dimensions. 
Nevertheless, Eq. (0) is still a crude approximation, because, generally speaking, more 
than one anomalous dimensions are expected to occur both for the twist-4 and the twist- 
6 terms. It is therefore clear that any phenomenological ansatz may introduce a model 
dependence in the final result, and in this respect it becomes important to try to estimate 
the theoretical uncertainties associated to the different higher-twist models as well as to the 
different perturbative approximations. Such a goal is beyond the aim of the present work 
and will be the subject of future investigations. 

In conclusion, we have analyzed the existing data on the (unpolarised) transverse 
structure function of the proton for large values of the Bjorken variable x by performing 
power correction analyses of the Nachtmann moments of Ref. || for 1 < Q 2 (GeV/c) 2 < 20. 
The results obtained adopting the NLO approximation and those including the effects of soft 
gluon resummation have been compared. We have shown that: i) higher twist extractions 
based on the NLO approximation for the leading twist are not reliable, and ii) the effects 
of the soft gluon resummation as well as the updated PDG value of a s [M\) should be taken 
into account at large x. 
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Figure 1. Ratio r n (Q 2 ) of the quark coefficient functions evaluated within the SGR ap- 
proach and at NLO (see Eq. (^)) versus Q 2 , adopting the values a s (M|) = 0.113 (a) 
and a s (M§) = 0.118 (b). Dot-dashed, dotted, short-dashed, long-dashed and solid lines 
correspond to n — 2, 4, 6, 8 and 10, respectively. 
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Figure 2. Nachtmann moments M n (Q 2 ) of the (unpolarised) transverse structure function 
of the proton versus Q 2 for n — 4 (a), n = 6 (b), n — 8 (c) and n = 10 (d). The open squares 
are the "experimental" moments as determined in ||, and the solid lines are the results of 
the twist analysis based on Eq. (HD adopting the SGR expression (|5]) for the leading twist. 
Long-dashed, short-dashed and dotted lines represent the twist-2, twist-4 and twist-6 terms 
separately, while the dot-dashed lines represent the total higher twist contribution. The 
error bars on the different twist contributions correspond to the uncertainties generated by 
the fitting procedure. The value a s (M|) = 0.118 is adopted. 
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Figure 3. The leading twist moments fi n (Q 2 ) resulting from the twist analyses of the 
Nachtmann moments of Ref. ||, based on the SGR expression (^) [solid lines] and on 
the NLO approximation (|j) [dashed lines]. Thin and thick lines correspond to the cases 
a s (Mf ) = 0.113 and a s (M|) = 0.118, respectively. The order of the moment is n = 4 (a), 
n — 6 (b), n = 8 (c) and n — 10 (d). 
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Figure 4. The same as in Fig. 3, but for the ratio of the total higher twist contribution 
HT n (Q 2 ) (Eq. (0)) to the leading twist term /i„(Q 2 ). 
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